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An inertial mass of a vortex can be calculated by driving it round in a circle with a steadily 
revolving pinning potential. We show that in the low frequency limit this gives precisely the same 
formula that was used by Baym and Chandler, but find that the result is not unique and depends on 
the force field used to cause the acceleration. We apply this method to the Gross-Pitaevskii model, 
and derive a simple formula for the vortex mass. We study both the long range and short range 
properties of the solution. We agree with earlier results that the non-zero compressibility leads to 
a divergent mass. From the short-range behavior of the solution we find that the mass is sensitive 
to the form of the pinning potential, and diverges logarithmically when the radius of this potential 
tends to zero. 



PACS numbers: PACS: 67.40.Vs, 47.37.+q 



INTRODUCTION 



Conflicting results on the mass of a quantized vortex 
in a neutral superfluid can be found in the literature. 
Popov [l[ and Duan Q argue that the mass per unit 
length is infinite, while Baym and Chandler [3] argue 
that it is negligible. We have developed a method for 
studying the problem of vortex mass and applied it to 
the Gross-Pitaevskii model [1, [j| for superfluids, which 
is known to describe weakly interacting bosonic atoms 
well, and is believed to give a useful qualitative descrip- 
tion of superfluid 4 He. We have obtained new analytical 
results and laid the basis for detailed numerical calcula- 
tions. Firstly, we show that our method is equivalent to 
the method described by Baym and Chandler. Secondly, 
we specialize to the Gross-Pitaevskii model and derive a 
new and compact formula for the vortex mass. Thirdly, 
we agree with Popov and Duan that there is a divergent 
contribution to the mass due to the expansion of the fluid 
at large distances from the vortex core. Fourthly, we get 
the new result that the mass is sensitive to the form of 
the pinning force that sustains the motion of the vortex 
in the presence of the Magnus force, and that the mass 
diverges as the range of this pinning force tends to zero. 
This suggests that an unambiguous vortex mass may not 
exist, and that inertial effects in vortex dynamics may be 
scenario-dependent Q. 

Vortex motion is important in many physical scenar- 
ios, from collective modes in neutron star matter to de- 
cay of superconducting currents via quantum nucleation 
of vortices and quantum turbulence in superfluids Q. 
The universal topological constraints on vortex behav- 
ior are an important asset in understanding these di- 
verse and important phenomena; vortex inertial mass is a 
salient uncertain parameter in the otherwise highly con- 
strained vortex dynamics. Since theories of a universal 



vortex mass based on different assumptions do not agree, 
our contribution here is to consider a specific scenario in 
which a vortex mass may be identified without ambiguity. 

To ascertain the mass of an object in the conventional 
way one determines the acceleration due to an applied 
force. Most of the methods that have been used so far 
to determine the mass of a vortex in a superfluid have, 
however, avoided the complexity of actual vortex accel- 
eration. In Popov's [l| work analogies between superflu- 
idity and electrodynamics are used, where vortices cor- 
respond to lines of electric charge and sound waves to 
electromagnetic waves. The energy E needed to form a 
vortex at rest is calculated, and the mass My is found 
by using the Einstein formula E = Myc 2 s , where c s is 
the sound velocity in the superfluid. The long range of 
the field due to a line charge leads to an infrared diver- 
gence of the self-energy, and so to an infrared divergence 
of the vortex mass. Duan gets similar results which he 
argues are due to the breaking of gauge invariance for 
the neutral superfluid. In the Baym and Chandler work 
[3j the calculation of mass is based on a calculation of 
the extra energy associated with a vortex forced to move 
with a constant velocity relative to the fluid. The authors 
argued that this energy is primarily associated with the 
core of the vortex, and that this is small, but they did 
not describe a detailed calculation. 

Accelerating vortices can be analyzed in some theo- 
ries, however. In classical incompressible hydrodynamics 
the mass of a hard-cored vortex can be calculated from 
the resonant frequency at which the vortex can move, 
without any applied force, on a circular orbit ll]. The 
mass obtained is just the finite mass of fluid displaced by 
the vortex core, plus the mass of the material contained 
in the core. In the Gross-Pitaevskii model, it has been 
shown that damping due to radiation of phonons is too 
strong to leave any trace of a similar resonance, unless 
the vortex core is loaded with a large additional mass 
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0, [l(| ■ We have carried out a more general and detailed 
calculation similar to that of Quist [10|, and we get sim- 
ilar results, even when the bulk modulus is much higher 
than it is in the Gross-Pitaevskii model Q. 

In the classical incompressible fluid, the same vortex 
mass can also be found by analyzing circular vortex mo- 
tion under an applied external force (acting on the core). 
A vortex moving with constant velocity v relative to the 
fluid must be acted on by a transverse force, equal to 
p s KV per unit length, where k is the circulation round 
the vortex and p s the superfluid density. This is the 
Magnus effect, which follows from conservation of mo- 
mentum and circulation, and so applies to vortices in 
quantum fluids as well as classical. We believe it may 
be an important omission in earlier work, that remarks 
such as "Now give the vortex an instantaneous velocity 
v" are made, without considering how the necessary 
force is to be applied. We therefore consider an external 
'pinning' potential, which repels the fluid particles, and 
hence attracts the low-density vortex core. By moving 
this potential around a circular orbit, we drag the vor- 
tex. We take it as our definition of the vortex mass, that 
the force required to do so will differ from the Magnus 
force by the centripetal force, equal to the product of 
vortex mass and acceleration. An orbiting vortex will 
generally excite waves, requiring an azimuthal force to 
maintain constant speed against radiative damping; we 
avoid this issue by considering slow enough orbits that 
this effect becomes negligible. 

We therefore begin by considering the perturbatively 
slow circular motion of a driven, pinned vortex in a ho- 
mogeneous quantum fluid, and by showing that our mass 
definition gives the Baym-Chandler formula |3J] in the low 
orbit frequency limit. 



PERTURBATION THEORY IN A ROTATING 
FRAME 



We consider a many-body hamiltonian H which in- 
cludes the kinetic energy of the atoms, the translation- 
invariant interaction Ui between them, and a moving pin- 
ning potential V p (r,t), which we take to have an instan- 
taneous axis of symmetry. We take it to move round 
a circle of radius a at a constant angular frequency lo. 
We ignore the confining boundaries that hold the sys- 
tem in place and hold the density constant. In a coordi- 
nate system that rotates with the pinning potential the 
Hamiltonian becomes time- independent; using polar co- 
ordinates centered on the axis of the pinning potential, 
rather than on the center of the circular orbit, we replace 
dt by dt+vd y +uid,p, where v — auj. We therefore consider 
an eigenstate of the rotating, shifted frame hamiltonian, 



where we have chosen units with h = 1. We write H' = 
Hr — ivdy and treat the term proportional to v as a 
perturbation. The vortex state for v = is l^o), with 
total energy Eq. To leading order in v we thus have 

\9 V ) = [1 + (Eo-HT 1 (ivd y )] \^ ) (2) 
E v = E + v 2 (y \dy(E - H'y'dy^o) ■ (3) 

The expectation value of the force exerted by the pinning 
potential is, to lowest order in v, 

F x = 2vS(V \(d x V p )(E - H')-%\9 ) ■ (4) 

We can use the commutation relation 

[d x ,H'] = d x V p + icjd y (5) 

to rewrite this as 



F x = iv(d x ip \dyipo) - it>(<9yVo|d x ^o) 
+2vu(ip \dy(Eo- H')-%\ip ) 



(0) 



The first two terms on the right can be shown, by apply- 
ing Stokes theorem, to give pv times the circulation [12J, 
and this is the standard form for the Magnus force. The 
coefficient of the acceleration vlj is our vortex mass My . 
Comparing it with (3) above, we recover the formula of 
Baym and Chandler [3j, 



d 2 E v 
dv 2 



(7) 



H + w(i<90 + 1) + ivdy 



(1) 



but now derived from the force on an accelerating vortex. 

Rather than proceeding further with the full many- 
body problem, we specialize to the Gross-Pitaevskii mean 
field theory 0, We study the asymptotic properties 
of the solution for large distances, and display the diver- 
gent contribution to the vortex mass for this compressible 
quantum fluid. We then study the short-range properties 
and consider the influence of the form of the pinning po- 
tential on the vortex mass. We find that the vortex mass 
depends strongly on the pinning potential, and diverges 
when its radius goes to zero. Finally, we briefly consider 
the implications of these results. 



GROSS-PITAEVSKII EQUATION AND 
COMPRESSIBILITY 



We consider the Gross-Pitaevskii equation in a frame 
of reference moving with speed v. which takes the form, 
in appropriately rescaled units [1J], 

- V 2 V + V p {r)i) + - 1)V> = -ivdyip . (8) 

The vortex solution for the limit v — > has the form 
studied by Gross and Pitaevskii 0, U|, ipo{ r ) exp («(/>), 
where i/'o satisfies 

- < ~ ;V>o + + (V p + \M 2 ~ l)^o = ; (9) 
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we take ipo to be real. We use the methods of Fetter [l5| , 
so that, to first order in perturbation theory, we write 
the solution as 



(10) 



m=0,2 



We can substitute this back into Eq. ([8]) to get the equa- 
tion, to first order in v, 



1 



-X' m + —Xm + {V p + 2il>l - 1 + i m cj) X r. 



+%x 



2-mJ 



im<t> 



■ < 



-ipo - ism4> ip' ; (11) 



we can take Xm to be real. In this paper we discuss the 
limit u) — ► of Eq. ©, and in this limit the Galilean 
invariance of Eq. © in the region where the pinning po- 
tential is zero gives a particular integral 



XQ = -\ r ^0 , XI = J^O 



(12) 



This satisfies the boundary conditions for small r, but a 
solution (/01/2) of the homogeneous equation 

-fm~ -fL + — fm + (V P + 2^ 2 - l)/ m + i>lh- m = 

(13) 

must be added to meet the boundary condition at r — > 00. 

For large r we use f± = /o ± fo, which satisfy the 
homogeneous equation 



= 



(14) 



where 



d 2 1 d 



++ 



dr 2 
d 2 



r 



dr 
1 d 



+ — + V p + 3^ 2 - 1 , (15) 



r dr 



The two solutions of the homogeneous equation bounded 
at large r fall off like /_ ~ 1/r and / + ~ exp(— y/2r)/y/r, 
and there is an unbounded solution with /_ ~ r, which 
can be used to cancel the dominant term in —ripo/2- We 
thus look for a bounded solution of Eq. (|TT|) of the form 



X+ = /+ , X- = f- - rV>o/2 



(16) 



To use the Baym-Chandler formula for the centrifugal 
force in Eq. ([5]), we can calculate the perturbed wave 
function to first order in v, and then calculate the excess 
energy that it contributes. This gives, with the use of 

Eqs. ©, mm, (USD, 



/+ /- - ^ripa 



#4 



/+ 



r dr , 



4tt 



' 1 f 

[(f---r^' + —Mrdr . (17) 



To find out whether the effects of the compressibility at 
large r give rise to a divergent vortex mass, as Popov [l[ 
and Duan Q have argued, we must find the asymptotic 
behavior of solutions of Eq. (Tl4]) for large r and substitute 
this into Eq. ([P?]). We have ip ~ 1 - l/2r 2 . Bounded 
solutions of the homogeneous equation have the form 

/--^ ~^ lnr+ S ; (18) 
the coefficient ai is determined by the small r bound- 
ary. The potential energy given by the /+ term gives 
an integrand proportional to r _1 , and so gives the loga- 
rithmically divergent expression for mass found by Popov 
[H and Duan 0. For real systems this integral will be 
cut off at large distances by the finite size of the system, 
the presence of other vortices, or by non-zero frequency 
effects. 



INFLUENCE OF THE PINNING POTENTIAL 

To study the influence of the pinning potential V p we 
take the specific case of a hard core repulsive potential 
of radius r c , so that ipo(r) = for r < r c . The particular 
integral given by Eq. (|12[) for the u> — ► limit vanishes 
at r = r c , so we must add solutions of the homogeneous 
equation that also vanish there in order to satisfy the 
boundary conditions for large r. To study these we go 
back to the representation used in Eqs. (fTTj) . (TT2")) . (fT3"j) . 
If r c << 1, two independent solutions for small r that 
vanish at r — r c can be constructed by combining solu- 
tions regular and irregular at the origin, and, close to r c 
where i/jq is small, these have the forms 



Jo 

f (2) 
JO 



-(ln^+ 7 )Jo(r)--Fo(r), /. 

0, A 2) *^Mr)+^Y 2 (r) 
ri 16 



(i) 







(19) 



Here J, Y are Bessel functions and 7 is the Euler constant. 
The solution to lowest order in i/}q is then 

ifr w e ^{^o+«[^osin0+a lX «e-^+a2X2 2) e # ]}, (20) 

where 01, 0,2 are real. 

Determination of these two coefficients requires a de- 
tailed integration of the homogeneous part of Eqs. (TT4")) 
and (|15p . but there is an important constraint placed on 
them by the fact that the Magnus force has magnitude 
2nv. For a hard core potential of radius r c the force per 
unit length can be written as 



F x =r< 



2tt 



\d r tp(r c , cj))\ 2 cos <j> d(j) 



(21) 
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From Eqs. ([20f and (fT9]) we get, at r = r c , 

RflK^e"**) =^o + — (ai + a 2 )cos^, (22) 

and the known value of the Magnus force gives 

f (r c )(a 1 +a 2 ) = l. (23) 

Energy considerations show that a 2 should be small, 
since a value of order unity would make a contribution 
of order r~ 2 to Eq. (fl7|) . so a\ must be close to ±1. The 
dominant contribution to Eq. (fl~7|) is 

M v « 4vr / ai lnr c J (r)(^o + — )dr , (24) 

where i? is a length of the order of unity. This shows that 
the vortex mass is sensitive to the form of the pinning 
potential, and diverges logarithmically as the core radius 
of the pinning potential tends to zero. It is the high 
quantum pressure near r c that leads to the large kinetic 
energy for small r c . 

We were actually led to this result by a different ap- 
proach in which we ignored the nonlinearity of the equa- 
tion for small r, and matched the solution to a region 
of incompressible flow for large r. This simplified model 
allowed us to get an exact expression of the vortex mass 
which displays this logarithmic divergence, and for which 
the wave function has the same general features at small 
r. This will be discussed elsewhere Q. 

DISCUSSION 

We have shown how to obtain an expression for the 
inertial mass of a stable quantized vortex in an infinite 
neutral superfluid by subjecting it to a straight, circu- 
larly symmetric, pinning potential V p which is slowly and 
steadily rotated about a parallel axis whose distance a 
from the vortex is large compared with the size of the 
vortex core. The vortex has a steady state in a frame 
of reference that rotates about this axis of rotation, with 
the same angular velocity u> that the pinning potential 
rotates. We use perturbation theory to study this steady 
state, and to find the force which the pinning potential 
exerts on the vortex to keep it in a steady rotation. The 
leading term in this expansion, proportional to v, gives 
the Magnus force, in a form which is closely analogous to 
the Magnus force acting on a moving vortex in classical 
fluid mechanics. The next term, proportional to vlu, has 



a coefficient that can be interpreted as an inertial mass 
of the vortex. 

It is essential to have a pinning potential to stabilize 
the position of the vortex in a rotating system, and it 
has to be strong enough to hold the vortex against the 
Magnus effect and the centrifugal force. We agree with 
Popov and Duan that the mass determined this way is 
logarithmically divergent in the low-frequency limit, and 
we have shown that it depends sensitively on the form 
of the pinning potential, diverging logarithmically as the 
radius of the pinning potential tends to zero. 

An important tentative conclusion of this work is that 
"the mass" of a vortex is not well defined, but depends 
on the process by which the mass is measured. We shall 
discuss this further in more detail elsewhere. 
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